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The study of light pseudoscalar quark-antiquark bound states in the Dyson-Schwinger approach 
with the effective QCD coupling enhanced by the interplay of the dimension 2 gluon condensate (A 2 ) 
and dimension 4 gluon condensate {F 2 ), is extended to the r\-rj complex. We include the effects 
of the gluon axial anomaly into the Dyson-Schwinger approach to mesons. The calculated masses, 
mixing and two-photon decay widths of y and rf mesons are in agreement with experiment. Also, in 
a model-independent way, we give the modification of the Gcll-Mann-Okubo and Schwinger nonet 
relations due to the interplay of the gluon anomaly and SU(3) flavor symmetry breaking. 
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1. INTRODUCTION 

The dimension-2 gluon condensate (A“ A a/ Q = (A 2 ) attracted the attention of some researchers well over a decade 
ago; e.g., see Refs, lj, y, J, 0, 0 ■ However, there was a wide-spread opinion that, since this condensate is not gauge 
invariant, it cannot have observable consequences and cannot play an important role in QCD. In contrast to that, 
the gauge-invariant, dimension-4 gluon condensate (F^F 0 ^) = (F 2 ) @| has, over 25 years now, been a subject of 
studies discussing its value (e.g., see Ref. Q) and implications for QCD. 

After it turned out more recently that the Landau-gauge value of (A 2 ) corresponds to a more general gauge-invariant 
quantity, it attracted a lot of theoretical attention §| 9, HU HU HU 113, HI, HU HU HU HU’ to quote just several of 
many papers offering evidence that {A 2 ) condensate may be important for the nonperturbative regime of Yang-Mills 
theories. In Ref. HUl: we argued that (A 2 ) condensate may be relevant for the Dyson-Schwinger (DS) approach to 
QCD. Namely, in order that this approach [19ll20l l2ll I22| leads to a successful hadronic phenomenology [which has so 
far been treated widely only in the rainbow-ladder approximation (RLA)], an enhancement of the effective quark-gluon 
interaction seems to be needed at intermediate spacelike momenta, Q 2 ~ 0.5 GeV 2 . Reference Q~5| showed that the 
interplay of the dimension-2 condensate (A 2 ) with the dimension-4 condensate (F 2 ) can provide such an enhancement. 
It also showed that the resulting effective strong running coupling leads to the sufficiently strong dynamical chiral 
symmetry breaking (D\SB) and successful phenomenology in the sector of light pseudoscalar mesons. In addition, 
the issues such as quark propagator solutions, p 2 -dependent dressed (“constituent”) quark masses and a more detailed 
discussion of the parameter dependence of the results were addressed in Refs. |23Ll24l . 

In the present paper, we extend the treatment of the nonzero-isospin light pseudoscalar mesons of Ref. E3 to 
the r\—rf complex. First, in the next section, the key result of Ref. [l^ . namely its gluon-condensate-enhanced 
interaction 03, is briefly re-derived in another, less rigorous and heuristic way. In the third section, we review how 
the DS approach employing such an interaction, can give the successful light (i.e., involving quark flavors q = u,d, s ) 
meson phenomenology, as this is also needed for the good description of the presently interesting 77 and rf mesons. 
Nevertheless, for rj and rf this is not enough because of the influence of the gluon anomaly and thus we explain how 
its effects are included in the manner of Refs. The implementation of the anomaly and the SU(3) flavor 

symmetry breaking, as well as their interplay, are presented independently of any concrete dynamics in subsection 
EU and in detail, because of some renewed interest in the flavor dependence of the mixing in the r\-rf complex (e.g., 
Ref. |27j] and references therein). After that, the masses, mixing angle and two-photon (77) decay widths of r\ and rf 
are calculated. Discussion and conclusion are in Sec. 0] 

2. STRONG COUPLING ENHANCED BY GLUON CONDENSATES 

Reference M showed how the interaction <E3, phenomenologically successful in DS studies in the Landau gauge and 
RLA, resulted from combining the form [Eq. 0 below] which the running coupling has in the Landa u-gauge DS studies 
[HilMliUlU an< ^ ideas on the possible relevance of the {A 2 ) gluon condensate [UIUIUIUI?ilU liH[lillTUITUIl4| . 
In the present paper, we give a simplified and more intuitive derivation thereof as follows. 
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The full gluon propagator D°^ v (k) in the Landau gauge is defined through the free gluon propagator D^ v (k )o and 


the gluon renormalization function Z[—kr) like this: 


D™{k) = Z{-k 2 )D^{k )o = Z(-k 2 ) — ( -g^ + 


kfj,ki, 


The full ghost propagator D G (k) is similarly defined by the ghost renormalization function G(—k 2 ): 

G(—k' 2 ) 


D G {k) = 


k 2 


The strong running coupling a s (Q 2 ) used in the Landau-gauge DS studies 0i III S! EE E| is defined 

a s (Q 2 ) = a s (g 2 )Z(Q 2 )G(Q 2 ) 2 , 


as 


(i) 


( 2 ) 


( 3 ) 


where a s (/i 2 ) = <? 2 /47t and Z(/i 2 )G(/i 2 ) 2 = 1 at the renormalization point Q 2 = /i 2 . Our convention is k 2 = — Q 2 < 0 
for spacelike momenta k. 

The functions Z and G can be expressed through the corresponding gluon (A) and ghost (G) polarization functions 
U A (Q 2 ) and n g (Q 2 ): 


Z(Q 2 ) 


1 


i , n A (Q 2 ) 
1 + Q 2 


G(Q 2 ) 


1 


i , n g (Q 2 ) ■ 
1 + Q 2 


( 4 ) 


Almost two decades ago, it was noted that in the operator product expansion (OPE) the gluon condensate {A 2 ) can 
contribute to QCD propagators; e.g., see Refs. [2, 3j|,@, 5j. Their (A 2 )-contributions to the OPE-improved gluon (A) 
and ghost (G) polarization functions were more recently confirmed by Kondo El For the Landau gauge adopted 
throughout this paper, three QCD colors (N c = 3) and four space-time dimensions (D = 4), their expressions for the 
polarizations become 


IL(Q 2 ) =m 2 + e>,;(l/Q 2 ), (i = A, G), 


( 5 ) 


m2 A = h = ’ (6) 

where via and m G are, respectively, dynamically generated effective gluon and ghost mass. Reference El estimated 
in A = 0.845 GeV (and found that it was, phenomenologically, a remarkably successful initial estimate) by using in 
Eq. © the lattice result q 2 (A 2 ) = 2.76 GeV 2 of Ref. a value compatible with the bound resulting from the 
discussions of Gubarev et al. |9j, 10] on the physical meaning of (A 2 ) and itspossible importance for confinement. 

As for the contributions Oi(l/Q 2 ) (i = A, G) in Eq. ©, it turned out @, 40 that they contain many kinds of 
mostly unknown condensates [e.g., gauge-dependent gluon, ghost and mixed ones, where terms oc (1 /Q 2 ) n [n > 1) 
were not considered at all]. The only practical approach at this point is therefore that these complicated contributions 
are approximated by the terms oc 1/Q 2 and parameterized, i.e., 

O a (1/Q 2 ) « ^ , O g (1/Q 2 ) « . (7) 

Thus, Ca and C G would in principle be free parameters to be fixed by phenomenology. However, as noted in Ref. 
EL in the effective gluon propagator proposed by Lavelle the Oa (1/Q 2 ) polarization is (for the Landau gauge 
and D = 4) given by the dimension-4 gluon condensate {F 2 ) as 


34 V c 7 ra s (F 2 ) (0.640 GeV ) 4 
~ 9(V 2 - 1)Q 2 “ Q 2 


Since Lavelle’s @2j propagator misses some unknown three- and four-gluon contributions @@ and since the precise 
value of a s {F 2 ) is still not certain, we regard the value Ca = (0.640 GeV ) 4 just as an inspired initial estimate. Still, 
together with the assumption C G = Ca , it was a very useful starting guess in our Refs. leading to very 

good phenomenological fits. 

We are now prepared to give a general, although heuristic argument why the contribution © of the dimension- 
2 (A 2 ) condensate to the gluon and ghost polarization functions ©, should indeed lead to the form of a e ff(Q 2 ) 
already found through a more detailed argument in Ref. fl5l| . Our first step is to assume that in the gluon and 
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ghost polarizations n A and 11(5, one can disentangle the perturbative (pert) from nonperturbative (Npert) parts, 
IR = II? ert + nf pert (i = A, G). At least for high momenta Q 2 , it is then possible to approximately factor away the 
perturbative from nonperturbative contributions; for i = A , 


Z{Q 2 ) 


1 + 


,pert 




(Q 2 ) 


1 + 


,Npert 


(Q 2 


= Z pert (Q 2 ) Z Npert (Q 2 ), 


( 9 ) 


where the approximation means neglecting the contribution of the term 

n pert (Q 2 )n* pert (Q 2 ) 

Q 4 


Analogously, 


G(Q 2 ) 


1 + 


T pert 


(Q 2 


1 + 


T Npert 


~w 


(Q 2 


= G pert (Q 2 )G Npert (Q 2 ). 


( 10 ) 


( 11 ) 


Since the general QCD coupling a s (Q 2 ) must reduce to the perturbative QCD coupling a per t(<3 2 ) for so very high 
Q 2 that nonperturbative contributions are negligible, Eq. 0 implies that 

a s (g 2 ) Z pert (Q 2 ) G pert (Q 2 ) 2 = a per t(Q 2 ) • (12) 

We can also assume, for high Q 2 , that nonperturbative parts are given by the OPE-based results of Refs. iiiiini 
In our present case, they amount to Eqs. (0-® for the gluon and ghost polarizations, and to the parameterization 
0 . This then gives 


Z Npert (Q 2 ) 


G Np ert(Q 2 ) 


1 < Q2 ~ r>4 


Q 4 


+Ca 

1 TF ' Q* 


(13) 

(14) 


Equations ®, ®, mil . fHJii . and mi) . considered together, then suggest an effective coupling a e g(Q 2 ) of the form 

a eS (Q 2 ) = otpert(Q 2 ) Z Npeit {Q 2 ) G Npert (Q 2 ) 2 . (15) 

Obviously, the above derivation of the coupling m is only heuristic, but we have already presented its more 
rigorous derivation in Ref. 0- 

In Refs. ®E!|2|, we discussed why, how and when the form 113 was sufficiently enhanced at intermediate Q 2 
to lead to the successful pion and kaon phenomenology when used in the DS approach in RLA. 


3. DS APPROACH AND ITS EXTENSION TO rj-p' COMPLEX 

In the DS approach to QCD, one solves the gap equation, i.e., the DS equation for quark two-point functions, 
namely dressed quark propagators 


S q (p) = 


1 

p/A q {p 2 ) - B q {p 2 ) 


Mp 2 )- 1 

yt- M q {p 2 ) 


(16) 


of various flavors q , and so explicitly constructs constituent quarks characterized by the dynamical masses A4 q (p 2 ) = 
B q (p 2 )/A q {p 2 ). The constituent quarks and antiquarks of respective flavors q and q' build meson bound states, which 
are solutions of the Bethe-Salpeter (BS) equation for the bound-state vertex T qq r. 

[T q9 >(k,P)]ef = J ^{S q ^+j)r qq '(^P)S ql (e-j)] gh [K(k-£)] h e f , (17) 

where e, /, g , h schematically represent spinor, color and flavor indices. Solving Eq. o for T qq ’ also yields M qq ’ , the 
mass eigenvalue of the qq' bound state. Unfortunately, the full interaction kernel K(k,£,P) for the BS equation is 
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TABLE I: Masses Mp, decay constants fp and 77 decay amplitudes Tp 7 of the pseudoscalar qq' bound states P = n,K and ss 
resulting from our gluon condensate-induced a e ff(Q 2 ) 1151 . along with the parameter values fixed in Ref. m by fitting the pion 
and kaon properties - see Eq. and the text below it. These masses and decay constants are the input for the description 
of the r\-r[ complex. The last column is the constituent quark mass A4,j(0) pertinent to the corresponding qq meson, i.e., 
Mu{0) = Mdifi) for P = 7 T and _M S (0) for P = ss. [Later, in Eq. 1401 . we will name the unphysical ss pseudoscalar meson r/g, 
but note that the mass M, ls 15111 , introduced in the NS-S mass matrix m . includes the contribution from the gluon anomaly, 
whereas M sq does not.] 


p 

Mp 

fp 

1/(4^ 2 /p) 

Tp 7 

M q ( 0) 

7V 

0.1350 

0.09293 

0.2726 

0.2560 

0.3842 

K 

0.4949 

0.1115 




SS 

0.7221 

0.1329 

0.1905 

0.08599 

0.5922 


not known. Also, the full kernels of the gap equations [which supply the quark propagators S q (l) (ITTH) to Eq. £□] 
involve the full gluon propagator and the full quark-gluon vertex, which satisfy their own DS equations. They in 
turn involve higher n-point functions and their DS equations, etc. This infinite tower of the integral DS equations 
must be truncated to make the problem tractable 1H| (20j, 12 ll 12211 . The approximations employed in the gap equation 
and the BS equation must be mutually consistent in order to preserve the important characteristics of the full 
theory. In the low-energy sector of QCD, the nonperturbative phenomenon of DySB is the most important feature. 
Phenomenological DS studies have therefore mostly been relying on the consistently used RLA, where DySB is well 
understood ,20. OS- 13 3|- The consistent RLA implies that for interactions between quarks one uses Ansdtze of 
the form 

[K{k)t 9 f = i47ra eff (-fc 2 Xjfc)o [y [y Y]hf (18) 

both in the BS equation 113 and the gap equation d for the dressed quark propagators ca, 

Sq\p) =rf-m q - 1471-/ (4^1 M-(p-Q 2 l^(p-Qoy7'‘S ? Wy 7" • (19) 

Here, fh q is the bare mass of the quark flavor q breaking the chiral symmetry explicitly. The case m q = 0 corresponds 
to the chiral limit where the current quark mass m q = 0. 

One solves the DS equation m for dressed propagators of the light quarks (q = u, d, s) using a, hopefully, 
phenomenologically successful interaction [in the present paper, it is given by Eq. G31- These light-quark propagators 
are then employed in BS equations for quark-antiquark (qq) relativistic bound states. 

The most important advantage of adopting RLA in the context of low-energy QCD is that D\SB and, consequently, 
the appearance of light pseudoscalar mesons as (almost-)Goldstone bosons in (the vicinity of) the chiral limit, is well- 
understood and under control in this approximation scheme nans 134 [ 35 ] ]. Solving the BS equation in the 
chiral limit gives the vanishing pion mass, M T = 0. More generally and precisely, the light pseudoscalar masses 
M qq oc sjrhqi as required by QCD through Gell-Mann-Oakes-Renner (GMOR) relation. All this is a manifestation 
of the correct chiral QCD behavior in the DS approach, in which all light pseudoscalar mesons (7r° , K 0± , K°, rf) 

manifest themselves both as qq bound states and (almost-)Goldstone bosons of dynamically broken chiral symmetry. 
This resolution of the dichotomy u qq bound state vs. Goldstone boson”, enables one to work with the mesons as 
explicit qq bound states while reproducing (even analytically in the chiral limit) the famous results of the axial 
anomaly for the light pseudoscalar mesons, most notably the n° — > 77 decay amplitude T 7 J = l/(47r 2 /^) in the chiral 
limit [ 33 ,[37]. This is unique among the bound state approaches - e.g., see Refs. (H| and references therein for 
discussion. Nevertheless, one keeps the advantage of bound state approaches that from the qq substructure one can 
calculate many important quantities (such as the pion decay constant f n ) which are just parameters in most of other 
chiral approaches to the light-quark sector. 

In the DS approach, the DySB is obtained in an, essentially, Nambu-Jona-Lasinio fashion, but the DS model 
interaction, encoded in various popular forms of a e g(p 2 ), is usually less schematic. Typically, as in Refs. mm 
13311391 i40l |4J for example, it combines the modeled nonperturbative component (which is strong enough to obtain 
required DySB) and the known perturbative QCD contribution. Our gluon condensate-induced a e g(Q 2 ) (1151) also 
behaves nonperturbatively at low Q 2 and perturbatively at very high Q 2 , as seen both in its earlier derivation [lfj 
and the alternative derivation in Sec. 0 Thus, as shown in Ref. [ljj, our interaction 115[l leads to the deep Euclidean 
behaviors of quark propagators consistent with the asymptotic freedom of QCD [ 33 . However, what enables the 
successful description of pions and kaons in Ref. fl^t and what is also crucial in the present paper, are the contributions 
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to the gap and BS equations at low and intermediate momenta. The crucial result is the behavior of the gap equation 
solutions for the mass functions M q (Q 2 ) at Q 2 = 0 to Q 2 ~ (500 MeV) 2 , where M q (Q 2 ) (q = u, d , s), due to DySB, 
have values consistent with typical values of the constituent mass parameters in constituent quark models. [When 
we need to be specific, we can, for definiteness, choose to call the Q 2 = 0 value A4 g (0) the constituent mass.] In the 
chiral limit, where, as usual in the consistent DS approach 19L i20t l2ll 1221 . we get the correct chiral Goldstone-pion 
behavior, our a e g(Q 2 ) 11511 with the parameters from Ref. Il5l . namely 

C A = (0.6060 GeV ) 4 = C G , m A = 0.8402 GeV, (20) 

gives us A4 u ,d(0) = 369 MeV, and empirically already quite acceptable values for the pion decay constant and the 
(qq) condensate. Away from the chiral limit (but still for isosymmetric u- and d-quarks), the same interaction yields 
A4 u ,d( 0) = 384 MeV [just 4% above Ai u ,d(0) in the chiral limit] with the bare mass rh Uy d = 3.046 MeV, for which value 
our model reproduces the experimental values of n 0 mass, f n + decay constant as well as n° —> 77 decay amplitude, and 
respects GMOR relation jTj|- The empirical values of the kaon mass and decay constant were also reproduced very 
well when one in addition takes fh K = 67.70 MeV [l5j]. The s quark constituent, dynamical mass is then A4 S (0) = 592 
MeV. (In our earlier papers [la, ya UM we also found that the results for m Ut d and fh s were rather robust. For 
example, the values quoted here because they are preferred when the interaction m is employed in the gap and BS 
equations, 11191) and CE3- are quite close to the values of m Uj d and m s preferred when the Jain-Munczek effective 
interaction [33J is used instead.) 

Up to accounting for the gluon anomaly, the results in Table [I] are the input from the well-described pion and kaon 
sector M which enables, without any refitting of the model parameters m u ,d, rn s and (HU), the good description of 
the 77 — 77 ' complex. 


A. Masses in the 7 — rj complex 

The description j25j, [2(| |42] of r\ and vf is especially noteworthy, as it is successful in spite of the limitations of 
the DS approach in the ladder approximation. For this description, the crucial issues are the meson mixing and 
construction of physical meson states. For the DS approach, they are formulated in Refs. (2f| j2(|, where solving of 
appropriate BS equations yields the eigenvalues of the squared masses, M 2 S , M?,, M 2 S , and M 2 S , of the respective 
quark-antiquark bound states | uu), \dd), |ss), and | us). The last one is simply the kaon, and M U g is its mass Mr-. 
Nevertheless, the first three do not correspond to any physical pseudoscalar mesons. Thus, M 2 s ,M 2 j, M 2 g do not 
automatically represent any physical masses, although the mass matrix (to be precise, its non-anomalous part , which 
vanishes in the chiral limit) is simply 


Mna ~ 


Muu 0 0 

0 M dd 0 

0 0 M 2 S 


( 21 ) 


in the basis \qq),(q = u,d,s). However, the flavor SU(3) quark model, and especially the almost exact isospin 
symmetry, leads one to recouple these states into the familiar SU(3) octet-singlet basis of the zero-charge subspace of 
the light unflavored pseudoscalar mesons of well-defined isospin quantum numbers: 


k°> = 

-^={\uu) - | dd)) , 

(22) 

\vs) = 

-^=(|7m) + | dd) - 2 ss)) , 

(23) 

1 m) = 

-^=(|7m) + | dd) + |ss)) . 

(24) 


With | uu) and \dd) being practically chiral states as opposed to a significantly heavier |ss), Eqs. lt2Hi do not define 
the octet and singlet states of the exact SU(3) flavor symmetry, but the effective octet and singlet states. However, 
in spite of this flavor symmetry breaking by the s quark, these equations still implicitly assume nonet symmetry in 
the sense ||i| that the same states | qq} (q = u, d, s ) appear in both the octet member r/ s lt!231) and the singlet ?? 0 IMJ- 
Nevertheless, in order to avoid the Ua(1) problem, this symmetry must ultimately be broken by gluon anomaly at 
least at the level of the masses of pseudoscalar mesons. 
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In the basis E2J> J2H . the non-anomalous part of the mass matrix of 7r° and etas is 


Mna ~ 


Ml 0 0 

0 M| g M| 0 
0 Mq 8 Mqq 


(25) 


showing that the isospin / = 1 state 7 r° does not mix with the 1 = 0 states 77 s and 770 , thanks to our working in the 
isospin limit, where M u n = M dd , which we then can strictly identify with our model pion mass M n . Since in this 
model we can also calculate Ml = (ss|M| rA |ss), this gives us our calculated entries in the mass matrix: 

Mi 8 = (Vs\M 2 na \vs) = \ {M% + \mI) = (594.72 MeV ) 2 , (26) 


Mi 0 = (vs\M 2 na \ Vo ) = M 0 2 8 = ^ (M 2 - M 2 S ) = -(487.05 MeV ) 2 


(27) 


Mio = (vo\Mna\vo) = l (\m 2 sS + Ml) = (431.22 MeV ) 2 . (28) 

The values on the far right of these equations were calculated from M v = M ua and M sS from Table |I] i.e., they result 
from the parameters fixed in Ref. m. These contributions are substantial, but if we take the chiral limit, all of them 
would tend to zero according to the GMOR relation, as required by the chiral symmetry of QCD. Thanks to this 
relation, even for the realistically large strange mass our approach has 

M% « 2 M 2 K - Ml (29) 

in a reasonably good approximation, at the 10 % level, whereby 

M^ » \m\ ~ \Ml ■ (30) 

Equation (EH) is thus revealed as a variant of the standard 77 s Gell-Mann-Okubo relation (® featuring only pion and 
kaon masses, which are not affected by the anomaly 1 . This is not surprising, as the role of M| g in Eq. E3 is the 
non-anomalous 77 s “mass”. Similarly, even Moo @, which by Eq. (El gives the “770 Gell-Mann-Okubo relation” 


Mio 


2.,2 1 


-Ml + -Ml. 
3 K 3 w 


(31) 


is just the non-anomalous part of the 770 “mass” M rj0 . It however requires the anomalous, chiral-limit-nonvanishing 
part to avoid the Ua( 1 ) problem. 

In our DS approach, all the model masses M q g> (< 7 , q' = u , d, s) and corresponding qq' bound state amplitudes are 
obtained in the ladder approximation. Thus, regardless of any concrete model, they are obtained with an interaction 
kernel which cannot possibly capture the effects of gluon anomaly. Fortunately, the large N c expansion indicates 
that the leading approximation in that expansion describes the bulk of main features of QCD. The gluon anomaly is 
suppressed as 1/N C and is viewed as a perturbation in the large N c expansion. It is thus a meaningful approximation 
|23 t° consider the gluon anomaly effect only at the level of mass shifts and neglect its effects on the bound-state 
solutions. 

In the chiral limit and, as it will turn out, the SU(3) flavor limit, the gluon anomaly is coupled only to the 
singlet combination 770 E3J- Only the 770 mass receives, from the gluon anomaly, a contribution which, unlike quasi- 
Goldstone masses M qq i's comprising . does not vanish in the chiral limit. As discussed in detail in Ref. |25j . 
in the present bound-state context it is most convenient to adopt the standard way (see, e.g., Refs. 0, 51]) to 
parameterize the anomaly effect. We thus break the Ua( 1 ) symmetry by shifting the 770 (squared) mass by an amount 


1 The kaon is protected from mixing not only by isospin, but also by strangeness. 
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denoted by 3/3 = X v (in the respective notations of Refs. [2(| [l(| and Ref. [25|). The complete mass matrix is then 
M 2 = M% A + M A , where 


M\ = 


0 0 0 
0 0 0 
0 0 3/3 


(32) 


The value of the anomalous 770 mass shift 3/3 is related to the topological susceptibility of the vacuum, but in the 
present approach must be treated as a parameter to be determined outside of our bound-state model, i.e., fixed by 
phenomenology or taken from the lattice calculations such as Refs. 0, E El !5J ■ 

We now want to incorporate the effects of the realistic breaking of the SU(3) flavor symmetry into the description 
of the gluon anomaly. At this point it is customary to go straight to the nonstrange-strange ( NS-S) basis H39ll - (I40L 
but before doing this, it is instructive to rewrite for a moment the matrix El in the flavor, \qq) basis, where it has 
the pairing form, 


Mi = /3 


1 1 1 
1 1 1 
1 1 1 


(33) 


since this may be the most transparent place to introduce the effect of flavor symmetry breaking into the anomalous 
mass shift. Namely, Eq. El tells us that due to the gluon anomaly, there are transitions | qq) —> | q'q'); q , q' = u, d, s. 
However, the amplitudes for the transition from, and into, light uu and dd pairs are expected to be different, namely 
larger, than those for the significantly more massive ss. To allow for the effects of the breaking of the SU(3) flavor 
symmetry, we can write 


(w\M\\q'c?) = b q b q > 


(34) 


where b q = a// 3 for q = u,d and b q = Xyffi for q = s. The anomalous mass matrix El is, in the flavor-broken case, 
thereby modified to 


(35) 


where the arrow denotes rewriting M\ in the octet-singlet basis El El - Comparison with Eqs. El and El shows 
that incorporating into the anomaly the flavor symmetry breaking, X 1, leads to the following. First, the anomaly 
influences the r] S <-> 770 transitions, reducing the negative Mf 0 El by /3 q^(2 — A' — A' 2 ). More notably, the rjg and 
r/o masses including both non-anomalous and anomalous contributions are given by Eqs. El, El, and El as 



1 

1 

A 


1 

O 

O 

0 

^ to 

II 

1 

1 

X 

-0 

0 1(1-X) 2 

^2(2 — X — X 2 ) 


A 

A 

X 2 


0 ^(2-X-X 2 ) b(2 + X) 2 


M 2 s =M 2 8 + -(l-X) 2 f3, 


(36) 


K= M oo + l(l + X) 2 P- ( 37 ) 

Not only M Vo is modified, but the interplay of the gluon anomaly and flavor breaking modified the Gell-Mann-Okubo 
relation El as the anomaly becomes coupled also to rjs and influences its mass M Vg . 

The Schwinger nonet formula, derived from the condition that the trace and determinant of M 2 = M% A + M A be 
equal to those of the same matrix in the basis of mass eigenstates, here M 2 = diag(M 2 , M 2 , M 2 ,), now acquires the 
new term on the right-hand side: 

(4 M% - 3M 2 - M 2 )(3M 2 , + M 2 - 3 M 2 K ) - 8 - M 2 ) 

= —4(M| - Ml) 3/3 (1 - X 2 ) , (38) 

were we also used Eq. El- The usual Schwinger formula is known to be satisfied well for the vector and tensor 
nonets, but not for the pseudoscalar nonet (dill . Equation El reduces to the usual Schwinger pseudoscalar-meson 
relation for the limit of no anomaly, 3/3 —> 0, but also for just X —* 1, the limit of no influence of the flavor symmetry 










breaking on the anomalous mass shifts. Thus, introducing only the anomalous shift of the rjo mass still yields the 
usual Schwinger relation, as noted by Ref. :5l| in a different approach. 

The pion remains decoupled from the etas as long as one stays in the isospin limit; i.e., after one adds the anomalous 
contribution M\ (1351 to Eq. (1251) . one still can restrict oneself to 2 x 2 submatrix in the subspace of etas. The 
calculationally convenient basis for that subspace is the so-called NS-S basis: 


1 _ - 1 2 
I Vns) = + \dd)) =—\r] S ) + -\rj 0 ) , 

I vs) = I s ®) = ~\Z§I 7? 8 ) + Tfl 770 )' 


(39) 

(40) 


In this basis, the r\-rj mass matrix is 


M 2 

VNS 

M 2 

VsVNS 


ML + 2(5 

s/2(3X 

_> 

M 2 0 

M 2 

L VNsris 

M 2 

VS J 


V2/3X 

Mf s + pX 2 _ 


l 

O 


where the indicated diagonalization is given by the NS-S mixing relations 2 


I V) = cos cj) \ t]ns) - sin</> \r] S ) , \rf) = sin <j)\rj NS ) + cos</> \r] S ) 


(42) 


rotating vns,VS to the mass eigenstates 77 , 77 '. (In the last section we will use the effective-singlet-octet mixing angle 
6, defined by analogous mixing relations where rj^s —»■ Vs, Vs VOiP Q- It is related to the completely equivalent 
NS-S mixing angle <f> as 6 = <f> — arctan s/2 = </> — 54.74°.) 

Now the NS-S mass matrix m tells us that due to the gluon anomaly, there are transitions \rjffs) |? 7 s)- As in 
the argument above Eq. (TT51) . the amplitude for the transition from, and into, VNSi need not be the same as those 
for the more massive rjs . The role of the flavor-symmetry-breaking factor X is to allow for that possibility. There are 
arguments ( 2 fj, supported by phenomenology, that the transition suppression is estimated well by the nonstrange- 
to-strange ratio of respective quark constituent masses, Ai u and A4 S . Due to the Goldberger-Treiman relation, this 
ratio must be close |25|, [26j to the ratio of rj ns and rjs pseudoscalar decay constants fr\ NS = fn and f vs = f sS . In 
other words, we can estimate the flavor-symmetry-breaking suppression factor as X ss A4 U /A4 S or X « fn/ fun- (Yet 
another, but again closely related way of estimating A', is from the ratios of 77 amplitudes, as in Refs. |26L l4o| A In 
the present paper, we use X ss fn/fast because these decay constants are not only calculable in the DS approach, 
but also (in contrast to “constituent quark masses”) defined without any arbitrariness, and in the case of even 
experimentally measurable. Our present model result fn/fss = 0.6991 (see Table UJ is reasonably close to A exp « 0.78 
extracted phenomenologically I26i] from the empirical mass matrix m 2 xp featuring experimental pion and kaon masses, 
or, after diagonalization, experimental p and rj masses - see Eq. (ECU) below. 


B. Mixing angle and other results from r\—v/ mass matrix 


In this subsection, let us first see what hints we get from phenomenology. In our present notation, capital M a ’s 
denote the calculated, model pseudoscalar masses, whereas lowercase m a ’s denote the corresponding empirical masses. 

From our calculated, model mass matrix GD we form its empirical counterpart m by i) obvious substitutions 
M u a = M 7 r —» m, r and M sS — > m S s, and ii) by noting that m S g, the “empirical” mass of the unphysical ss pseudoscalar 
bound state, is given in terms of masses of physical particles as m 2 s ~ 2 m 2 K — m 2 due to GMOR. Then, 


to 2 + 2/3 s/2(5X 

_^ 

TO 2 0 

s[2(5X 2m 2 K — ?n 2 + (3X 2 

0exp 

1 

O 

^3 tO 

1 _ 


(43) 


where the arrow indicates the diagonalization (PI) for the angle value 4 > exp . 

Since M vM , obtained by solving the BS equation, is identical to our model pion mass M n , it was fitted to the 
empirical pion mass m x in Ref. |l5J. Similarly, M U g = Mk is fitted to the empirical kaon mass mpc, so that Eq. m 
implies M 2 5 ss 2 m 2 K — to 2 . We thus see that in our model mass matrix, the parts stemming from its non-anomalous 


2 The relation between the present approach and the two-mixing-angle scheme is clarified in the Appendix of Ref. J2jj]. 
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part M% a are already close to the corresponding parts in m 2 xp . We can thus expect a good overall description 
of the masses in and rf complex. We now proceed to verify this expectation. 

The anomalous entry /3, along with X (and then necessarily also (f>), is fixed phenomenologically if we require that 
they fit the masses m v and m v i in the empirical matrix (THU) to their experimental values. This is achieved by requiring 
that trace and determinant of irij Xp have their experimental values, which leads (e.g., see |53|) to the relations 


/^exp 


(m^ - ml)(m}, - ml) 
4 {m 2 K — m 2 ) 


= 0.2785 GeV 2 . 


(44) 


X K 



(to)), — 2 TO 


K + m l)( 2m K ~ m l - 


{nil, - nil){nil - nil) 


= 0.7791. 


(45) 


</> exp = arctan 


( to 2 , - 2 m 2 K + ml){ml 
{2m 2 K - ml - to;)) (to 2 , 



= 41.88°. 


(46) 


Now we want the analogous results from our theoretical mass matrix G2>. where only j3 is not a calculated quantity. 
For example, in this subsection we use X = f^/fss ~ 0.6991 from Tabled We thus fix only /3 by requiring that the 
experimental value of the trace ?u 2 + to 2 , ~ 1.22 GeV 2 be fitted by the theoretical mass matrix SB- This yields 


P = 2+X 2 t + ^ XP “ + M ] ’ ( 47 ) 

whereby X, M n o, and M sS from Tabled give us j3 = 0.2723 GeV 2 , in good agreement with /3 exp (EH) obtained from 
the empirical mass matrix m- 

Since M 2 n + M 2 g = 2 m 2 K holds to a very good approximation due to GMOR, our approach satisfies well the first 
equality [from Eq. 1471) ] in 


P (2 + A' 2 ) = to 2 + to 2 , - 2 m 2 K = x , (48) 

where the second equality is the Witten-Veneziano (WV) formula Q, with y being the topological susceptibility of 
the pure Yang-Mills gauge theory. The WV formula with experimental masses and jl yields % ss (179 MeV) 4 , which 
is in good agreement with our value x = (177 MeV) 4 , implied by Eq. iTISl) and our model values of X = f„/f^ and 
/3 from Eq. G3- Our prediction for \ is a l so iu reasonable agreement with the recent lattice results j47l l48l I49l | 
considered in Ref. ^3- The central value of Lucini et al. Xi = (177 ± 7 MeV) 4 0 agrees precisely with our value, 
while their value obtained by a different method, X 2 = (184 ± 7 MeV) 4 , is higher but still consistent with our y. 
However, = (191 ± 5 MeV) 4 of Ref. 0| is too high for that. On the other hand, our result is only marginally too 
hig h to be consistent with the most precise lattice topological susceptibility so far, Xi = (174.3 ± 0.5 ± 1.2lo'.2 MeV) 4 
l49l . In summary, our x is consistent with 

X = (175.7 ± 1.5 MeV) 4 , (49) 

the weighted average of the recent lattice results H3, |49|. Let us pause briefly to note that (3 does not have to be 
treated as the parameter to be fixed by fitting the masses, since Eq. GEJ> enables one to determine (3 from the lattice 
results on the topological susceptibility. In fact, this is how we get the results in the column C of Tabic ITT1 in the last, 
concluding section. 

However, now we continue with (3 from Ecp m fitting (to^ ~\~ m v ,) exp . For this value, {3 — 0.2723 GeV , the values 
of our calculated t)ns and Vs masses are 

M 2 = M 2 s + 2 (3 = M 2 + 2/3 = 0.5628 GeV 2 = (750.2 MeV) 2 (50) 


and 


M 2 s = M% + (3X 2 = 0.6545 GeV 2 = (809.0 MeV) 2 . (51) 

These are reasonable values, in good agreement with, e.g., tjns and ijs masses calculated in the dynamical SU(3) 
linear a model j46j . 
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The mixing angle is then determined to be <j> = 40.17° (or equivalently, 0 = —14.57°), for example through the 
convenient relation 


tan 2 cf> 


2y/2px 

M 2 - M' 2 

Vs lv± VNS 


The diagonalization of the NS-S mass matrix gives us the 77 and 7/ masses: 

M 2 = cos 2 ^ M 2 NS - V2/3X sm2(l) + sin 2 <P M 2 S 
M 2 , = sin 2 cj> M 2 ns + V2(3X sin 2 <f> + cos 2 <j> M 2 s 


(52) 


(53) 

(54) 


Plugging in the above predictions for /3, X , M rjNS . and M vs , our model 77 and r( masses then turn out to be M v = 579 
MeV and M v i = 939 MeV, in good agreement with the respective empirical values of 547 MeV and 958 MeV. 

The model values predicted in this and in the next subsection are summarized in column A of Table un Columns 
B and C give analogously obtained results, but with differently chosen either flavor breaking parameter X or the 
topologically susceptibility, i.e., /3. The comparison of our predictions shows they are robust under these variations. 


C. Two-photon decays of and r( 

Having obtained the predictions for the mixing in the 77-77' complex, we can get the predictions for the 77 decays 
of 77 and rf from the decay amplitudes T 7 0 7 and JT 7 already given in Table [I] However, for the sake of completeness, 
let us first briefly review how these amplitudes are obtained. 

The transition between the neutral pseudoscalar meson P and two photons 7 (fc) and 7 (k!) with momenta k and k! 
can be described by a scalar amplitude we denote Tp(fc 2 ,/c' 2 ) [3£j. The special case of the decay P —> 77 into two 
real, on-sliell photons corresponds to the k 2 = k' 2 = 0 amplitude 

T 77 = T P (0,0) = const, (55) 

so that by integrating over the phase space and summing over the photon polarizations one gets the decay width 

r(P- 77 ) = ^^|T 77 | 2 , (P = 7T°, V ,r,'). (56) 

The calculation of the electromagnetic transition amplitudes proceeds in the same way as in our earlier papers 
such as Refs, [lj, |25|, [2(| yij HQk EH since the incorporation of the quark-photon interactions is the same 

as adopted there through the scheme of a generalized impulse approximation, where all propagators, bound-state 
vertices, and quark-photon vertices are dressed. (In the present application, this impulse approximation is illustrated 
by the qq pseudoscalar -77 triangle graph in Fig. ^) They are all dressed mutually consistently, so that the pertinent 
Ward-Takahashi identities are respected (e.g., see Refs. ). This is necessary for reproducing exactly and 

analytically anomalous 77 (on-shell) amplitudes 3 in the chiral limit, and requires the usage of a dressed quark-photon 
vertex satisfying the vector Ward -Takahashi identity. We employ the Ball-Chiu vertex [6(j, which is widely used 
(e.g., see Refs. [2(| El El) EE EH EE [U EH and references therein). 

If one works in the NS~S basis m - m , the diagonalization (B 21 i and the amplitudes 

TXs = + = jT7, (57) 

77 = TV, (58) 

give the (mixing-dependent) amplitudes of the physical particles 77 and 77': 

T 77 = cos ct> T 77 s - sin <j> T 7 S 7 , (59) 

T 77 = sin 0 T 7 J S +cos0 T 77 . (60) 


3 And others, notably the “box anomaly” process 77r 


Tr+Tr 0 ; see Refs. 

mmEnn 
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FIG. 1: The triangle graph (and its crossed graph) relevant for the interaction of two photons of momenta k and k' with 
the neutral pseudoscalar meson P of momentum p represented by the quark-antiquark bound-state vertex Tp(q,p). The 
quark-photon coupling is in general given by dressed vector vertices r M (qi, $2), which in the free limit reduce to where 

Q = diag(Q„, Qd, Q s ) = ediag(+2/3, —1/3, —1/3) is the flavor SU(3) quark charge matrix and e is the electromagnetic charge, 
e = \J 47T fte m . 


In the preceding subsection, describing the breaking of the SU(3) flavor symmetry by X = f^/f s g led to cf> = 40.17°. 
This mixing angle value yields the physical amplitude values 

T 77 = 0.2706 GeV " 1 , T 77 = 0.3410 GeV" 1 , (61) 

which agree very well with the experimental amplitudes 

(T 77 ) exp = 0.2724 ± 0.0069 GeV -1 , (T 7 , 7 ) exp = 0.3417 ± 0.0060 GeV -1 . (62) 

The corresponding calculated and experimental widths T(P —► 77 ) (P = 77 , 7 /) are respectively displayed in the 
columns A and £ of Table QH As already mentioned, other columns give our predictions based on somewhat different 
choices for either the flavor symmetry breaking parameter A' or the topological susceptibility. Similarly to the results 
of the previous subsection, our results on 77 decays turn out to be robust under these variations, since the resulting 
changes in the mixing angle are not excessive. 

It is anyway very satisfying to note that the mixing-independent combination of our theoretical 77 decay amplitudes 

|T 77 | 2 + |T 77 | 2 = |T 7 ; S | 2 + |T 77 | 2 = 0.1895 GeV " 2 (63) 

matches very well the corresponding experimental value 

|(T 77 ) exp | 2 + |(T 77 ) exp | 2 = 0.1909 ± 0.0056 GeV " 2 , (64) 

because this match does not depend on the mixing angle at all. 


4. DISCUSSION OF RESULTS AND CONCLUSION 
A. On modification of analytic mass relations 

The interplay of the gluon anomaly and flavor symmetry breaking modifies some flavor SU(3) mass relations. The 
details of forming the mass matrix in subsection EH revealed the corrected rjs Gell-Mann-Okubo relation ESI) , its 
770 analogue m, and Schwinger formula (|38|) . It is obvious that such results are independent of the DS approach 
and even of any concrete dynamics, and depend only on the way the gluon anomaly and the SU(3) flavor symmetry 
breaking are implemented. It is thus not surprising that our Eqs. ESI), EZ1> , and ESI) are respectively equivalent to 
Eqs. (13), (15), and (38) of, e.g., Ref. J27] which does not address any dynamics. 

Our primary aim concerning these mass relations is to point out that all of them have been implicitly contained in 
our DS approach to the 77 - 77 ' complex and its results ever since Ref. |H] incorporated into the DS approach the effects 
of the flavor symmetry breaking on the anomalous mass shift. However, it is instructive to consider these relations 
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explicitly. Already Eq. 11351) showed that the usual Schwinger formula, which is badly violated for pseudoscalars, 
acquires the correction proportional to 3/3(1 — X 2 ). Note that the corrected Schwinger formula (13811 is identically 
satisfied if /3 and X are respectively given by Eqs. (1331) and 63), since these equations stem from the determinant 
and trace conditions, just like Eq. (I351i itself. As for the r]g Gell-Mann-Okubo formula, the original one (1331) was 
modified in Eq. (I331i by the term 


6 = 


|(1-A) 2 /3 = x 


ANr (1 - X) 2 

m 2+x2 


(65) 


where the last equality comes from expressing /3 by the topological susceptibility % through the WV relation (1331) . 
This result is independent of our concrete model, or even of the DS approach in general, as essentially the only 
strong simplifying assumption was the one of the (flavor-broken) nonet symmetry. If X = f^/ fss is chosen, the 
first-order flavor symmetry breaking estimate of f S g (which is in our model satisfied better than 0.5% - see Table 0, 
f S s ~ 2 Jk — fn , permits expressing the correction (El exclusively through experimental and lattice values. So, the 
experimental /^,/if and the central value of the weighted average l!351) of the recent lattice results UllEHE!] on X 
gives 5 ss (127 MeV) 2 independently of any model but in good agreement with our model result 5 = (128 MeV) 2 
[from our model values of X = f n /fss and /3 from Eq. (1371) ]. Equation JHHJ then yields M Vs = 608 MeV, in excellent 
agreement with the chiral perturbation theory result of some 610 MeV j6l|. It is satisfying that regardless of whether 
we take the lattice results for % and empirical results for decay constants, or our model values for /3 and X, our 
correction © is in reasonable agreement with Eq. (13.9) of Hagiwara et al. {62j (where our 5/M| s = A in their 
notation). It gives 6 = —14.1°, in good agreement with our model results discussed in the next, main subsection of 
the conclusion. 


B. Model results 

The bulk of our results on the 77 — 77 ' complex are summarized in Table GH namely our model predictions for the 
masses, the flavor symmetry breaking parameter A', then 3/3, the anomalous squared mass of 770 in the chiral or SU(3) 
flavor limit (where r\ = 770 ), the mixing angle 6 = (f> — arctan \/2 and the 77 , rj —> 77 decays. The results in the column 
A were also given in the text in the course of explaining our procedure. This column was obtained by assuming that 
the flavor symmetry breaking parameter X is given by our model value of f^/fss- The requirement (1471) that the 
inclusion of the anomalous mass contribution fit the sum of the (squared) experimental masses, (t? 7 2 + m 2 r ) exp , then 
fixes the mass matrix and the mixing angle which diagonalizes it. This angle is then used to calculate the 77 decay 
amplitudes and, equivalently, the corresponding widths. The values in the column B were obtained in the same way 
except that the flavor symmetry breaking parameter X is estimated from the ratios of our calculated 77 amplitudes, 
X = (TJ7/ Q 2 )/(Tgg / Q 2 ) , as explained in Ref. 26]. In the column C we again use X = f n /f sS as in the column A, 
but instead of fixing /3 by fitting (m 2 + m 2 ,) expi we obtain /3 through the WV relation (1331) from the central value of 
the weighted average 13*31) of the recent results on the topological susceptibility calculated on the lattice (43, E 3 0 ] ■ 
Thus, the column C gives the masses, mixing and 77 decay widths without any parameter fitting whatsoever. In 
spite of this, the results in the column C are just as consistent with the experimental masses and decay widths as the 
results in the columns in which /3 was used for fitting. 

In the column T> we give the best \ 2 ~fit to the experiment which our theoretical 77 decay amplitudes T^o and T 7 g 7 
can give using the mixing angle as a free fitting parameter, regardless of the results on the mixing angle from the 77-77' 
mass matrices. Nevertheless, the comparison with other columns in Table ITTI shows that what we find from 77, rj —> 77 
processes is actually close to what we find from the mass matrix, which is of course very satisfying. Actually, the 
comparison of all results in Table GH shows generally that our theoretical results in all columns are similar among 
themselves, exhibiting robustness under input variations, and all are in reasonable agreement with the experimental 
results in column £. 

Let us now compare the present results with our earlier work |25l (26 ] on the 77—77' complex. It also employed the 
consistently coupled DS approach but using the Jain-Munczek effective interaction 33 instead of our Eq. (1151) . The 
treatment of the 77-77' complex used here was largely formulated already in Ref. |25l , except that it did not consider 
the interplay of the SU(3) flavor symmetry breaking and the gluon anomaly, which was taken into account in Ref. |2G1| 
and found important for the successful meson phenomenology. The present paper, employing a different interaction, 
confirms its importance and further clarifies it by displaying explicitly the interplay of /3 and X in the mass relations 
(13611 - I3SI) . While the reproduction of the 77 and 77' masses in Ref. |26i| was similarly successful, we find that our 
77,77' —> 77 decay widths agree with the experiment much better now, although the error bars shrunk as the precision 
of the experimental widths substantially increased in the meantime {63j . 

Concerning the agreement with other approaches, we may point out that the 77-77' mass matrix obtained on lattice 
by UKQCD collaboration j64[ agrees reasonably well with our model mass matrix if we insert in it our values of /3, 
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TABLE II: Various theoretical results on 77 and rf and comparison with their experimental masses and 77 decay widths. All 
calculated quantities were obtained with parameters which gave the very good description of pions and kaons in Ref. El, 
i.e., Eq. 1201 and m u ^d = 3.046 MeV, m s = 67.70 MeV. The contents of the columns is this: A) Summarized predictions of 
subsecs. EH and l3Bl i.e., with X = fn/fss- B) X estimated from P —> 77 amplitudes (see text). In both A and B, the 
parameter (3 is fixed by fitting the experimental value of the mass matrix trace, Eq. 1471 . C) X = fn/fas again, while j3 
is not a fitting parameter but obtained from \ = (175.7 MeV) 4 , the weighted average 1401 of the recent lattice topological 
susceptibilities 1471 1481 1491 . The results of this column are thus obtained without any free parameters. T>) Starting from the 
calculated 77 amplitudes Tff and T ± 7 in Table [I] the mixing angle value is obtained as the fitting parameter in the best y 2 fit 
to the experimental amplitudes 1621 of 77 and 77 '. £) experimental values. The widths T(r/, rf —> 77 ) are calculated using the 
experimental masses. The 77 decay widths are in units of keV, masses are in units of MeV, 3/3 in units of MeV 2 , while X and 
the mixing angles are dimensionless. 



A 

B 

C 

V 

£ 

e 

-14.57° 

-14.93° 

-15.18° 

-14.27° 


M v 

579.3 

576.1 

577.1 


547.75 ±0.12 

M v , 

939.0 

941.0 

932.0 


957.78 ±0.14 

X 

0.6991 

0.7124 

0.6991 



3(3 

816951 

810835 

798060 



r( 7 ? -> 77) 

0.5034 

0.5114 

0.5170 

0.4968 

0.510 ±0.026 

r( 7 ?' -> 77) 

4.272 

4.229 

4.199 

4.308 

4.29 ±0.15 


X , M un , and M sS . We should also recall that already Ref. [2fil | clearly showed that our DS approach and results are 
not in conflict, but in fact agree very well with results in the two-mixing-angle scheme (reviewed and discussed in, 
e.g, Ref. :65|). Actually, Ref. ,26| showed that our results can also be given in the two-mixing-angle scheme, but it 
is defined with respect to the mixing of decay constants, and therefore in our case, as in the DS approach in general, 
the scheme with one angle defining the mixing of the states is more convenient. 

Another important general feature of the consistent DS approach which also holds in the present paper, is the 
correct chiral QCD behavior, so that all light pseudoscalar mesons constructed in our approach are both qq bound 
states and (almost-)Goldstone bosons of DySB. In the chiral limit for all quark flavors, the anomalous mass matrix 
M\ [with A' —7 1, Eq. (132ll ] is the only nonvanishing contribution as M^ A 0- Thus, in this limit 77 —7 rjg with 
vanishing mass, and rf — > 770 with the mass \/3j3 (= 0.904 GeV in column A and similarly in the other columns), which 
is the only non-vanishing pseudoscalar mass in that limit, being induced purely by gluon anomaly. This chiral-limit 
mass y/3j3 is thus only 6 % below the experimental rf mass. Even with the realistic breaking of the chiral symmetry 
(as in Table ITU which rises the 77 by more than 0.5 GeV, our mixing angles show that 77 and rf are much closer to, 
respectively, the (almost-)Goldstone octet 77 s and non-Goldstone singlet 770 than to tins and rjs (which is why in this 
section we switched from cf> to 9). For example, column A, with 6 = —14.57°, implies 

|r?) = 0.968 |tt 8 > + 0.257 |tto) , \rf) = - 0.257 |t») ± 0.968 It*,) , ( 66 ) 

where squaring of the coefficients cos0 = 0.968 and sin0 = —0.257 says that 77 is 93.4% r/% and 6 . 6 % 770 , and reverse 
for rf. 

In conclusion, the consistently coupled DS approach with the effective interaction O enhanced by the gluon 
condensates gives a very good description of the iq-iq 1 complex. This was achieved along the lines formulated in Refs. 
[25l tSfj , where the consistently coupled DS approach was extended by assuming the mass shift of the singlet 770 due to 
the gluon anomaly. After Ref. 0 found the model parameters for which the gluon-condensate-enhanced interaction 
dm leads to a sufficiently strong DySB, pions and kaons as (quasi-)Goldstone bosons of QCD, and their successful 
DS phenomenology, this minimal extension was the only new element in the otherwise fixed model. This was enough 
to successfully model 77 and if mesons without any parameter re-fitting, which is not very surprising after the success 
with pions and kaons in Ref. 0 - In Refs. [25| and especially [26 we have already given good descriptions of 77 and rf 
employing the Jain Munczek effective interaction 1.331 . which is however purely modeled at the low and intermediate 
energies. In contrast, the important intermediate-momentum behavior of the presently used interaction 113 may be 
actually understood in terms of gluon condensates, instead of just modeled. In addition, this interaction has presently 
given the description of 77 and rf which is on the whole somewhat better than in Refs, pfl l26l|. especially in view of 
the increased precision of the 77 , rf —> 77 measurements [63j. 
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